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I. INTRODUCTION 

The nucleon-nucleon (NN) interaction has been studied for more than 50 years, but 
it is still not fully understood. The research program proposed by Japanese physicists 
around 1950 [Q] and based on the idea that the outer part of the interaction is due to 
mesonic exchanges proved to be very successful. Quite generally, the spatial features of a 
given process are determined by the mass exchanged in the t-channel and lighter systems 
correspond to longer interaction ranges. The lightest NN exchange, associated with the 
one-pion exchange potential (OPEP) became a consensus in the sixties, giving rise to a 
generation of models in which waves with L > 5 were treated theoretically ||. 

The second layer of the interaction is much more complex and was studied in the following 
two decades, by means of a detailed treatment of the two-pion exchange potential (TPEP) 
||f|]. This process depends on an intermediate pion-nucleon (7rN) scattering amplitude and 
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reflects strongly its dynamical content. 

In the present decade, most of the research work on the NN potential was aimed at 
constructing the TPEP in the framework of chiral symmetry. This symmetry was devel- 
oped around 1960, for systems in which pions and nucleons were considered as elementary. 
Nowadays one believes QCD to be the basic theory of strong interactions and, accordingly, 
that pions and nucleons are made of light quarks, interacting by means of gluon exchanges. 
As the QCD Lagrangian predicts that gluons can interact among themselves, calculations 
at low and intermediate energies are very difficult. The usual strategy for overcoming this 
problem consists in working with effective theories that treat pions and nucleons as elemen- 
tary and include, as much as possible, the main features of the basic theory. The fact that 
the masses of the u and d quarks are close to each other and very small in the hadronic 
scale means that QCD is approximately invariant under the chiral group SU(2)xSU(2). One 
therefore requires the effective theories to possess approximate chiral symmetry, besides the 
usual Poincare invar iance. 

In low energy processes, chiral symmetry is realized in the Nambu-Goldstone mode and 
the vacuum is filled with a condensate, that allows the excitations of massless collective 
states, identified with the pions. The breaking of the symmetry, due to the quark masses at 
the fundamental level, is associated with the small pion masses in the effective theories. 

Chiral symmetry is very important in the theoretical treatment of two-pion exchange 
because it constrains the intermediate 7rN amplitude. At low and intermediate energies it 
is given by a nucleon pole contribution, superimposed to a smooth background 0. The 
symmetry is responsible for large cancellations within the nucleon sector that, at once, 
settle the scale of the problem and amplify the role of the background. The latter is very 
important, since the chiral nucleon sector in isolation does not suffice for explaining 7rN 
experimental data. 

The construction of the TPEP in the framework of chiral dynamics motivated most 
of the research of this decade, beginning with the work of Ordonez and van Kolck ||, who 
considered a system containing just pions and nucleons. Several works followed, dealing with 
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complementary aspects of the problem |7]|| and nowadays this part of the NN interaction 
is well understood. Predictions for NN observables produced by just the OPEP and the 
chiral TPEP, assumed to represent the full interaction for distances larger than 2 fm, were 
calculated and shown to agree well with experiment ||. Therefore the effort based on chiral 
symmetry led to an important refinement of outer part of this interaction and brought 
theoretical constraints to waves with L > 3. 

In the case of the NN interaction, it is important to note that the importance of chiral 
symmetry depends strongly on the process one is considering. In the case of the OPEP, 
for instance, it is completely irrelevant, for predictions from chiral 7rN Lagrangians coincide 
with those arising from interactions without any symmetry. All Lagrangians, symmetric 
and non-symmetric, produce exactly the same basic pion-nucleon vertex, showing that chiral 
symmetry is compatible with and, at the same time, irrelevant for the OPEP [|1(J. In the case 
of the TPEP, on the other hand, the symmetry is crucial. It produces internal cancellations 
in the intermediate 7rN amplitude which yield a potential that vanishes in the chiral limit. 

To our knowledge, only contributions due to the exchanges of one and two pions have 
been so far studied in the framework of chiral symmetry^]. In order to extend this picture, 
here we study the component of the NN potential due to the exchange of three uncorrelated 
pions. This system has a mass around 450 MeV and its effects should be longer than those 
of the vector mesons usually present in one-boson exchange potentials. The basic interaction 
is closely related to the amplitude for the process 7rN— > tnrN. Hence, we review briefly the 
main features of this reaction in sect. II and calculate the potential in sect. III. 



II. PION PRODUCTION 

The NN interaction mediated by the exchange of three uncorrelated pions is derived from 
T c ba, the amplitude for the process n a (k)N(p) — > 7i b (q) 7T c (q')N(p'). This amplitude is given 



1 For an early work on problem see, for instance, ref. [11 



by the diagrams of fig. 1 and written as the sum of T^ ba , representing the class of processes 
with a pion pole in the t-channel, and a remainder, denoted by T c f, a , whose dynamical 
content is indicated in fig.l. In the framework of chiral symmetry, this last amplitude is 
given by a basic family of diagrams, involving only pions and nucleons, supplemented by 
other processes, containing deltas, rhos, omega and the tcN a-term. 
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FIG. 1. Diagram for pion production. 
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There are many alternative ways of implementing chiral symmetry. In particular, the 
subset of diagrams in the pure pion-nucleon sector, corresponding to the minimal realization 
of chiral symmetry in this problem, may be evaluated by means of non-linear Lagrangians 
with 7rN couplings that may be either pseudovector (PV) os pseudoscalar (PS). Denoting 
the pion field by 0, and defining / = (f% — <fi 2 ) 1 / 2 , we have 



PV 



W+f*) 2 
C PS = jC« + Ni0N- gN(f + it ■ <f> l5 )N 



where 



■■(dpfPf+d^-dvfl+Uiff 
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In these expressions, ip and N are the nucleon fields with non-linear and linear transformation 
properties, /i and m are the pion and nucleon masses and /„., g and Qa are respectively the 



pion decay, the 7rN coupling and the axial decay constants. It is important to stress that 
these Lagrangians, in spite of their different aspects, have the same dynamical content and 
physical results do not depend on the particular version one adopts, as demonstrated on 
general grounds [[12]. 

The pion production amplitude has the general form 



iTrhn = U 



S bc r a {A* + A) + 5 ac r b {B* + B) + S ab r c (C7* + C)+ ie cba E 



75 u 



(4) 



where the tags 7r and bar refer to the pion-pole and background contributions. At threshold, 
this amplitude is usually written as 



l th 



cba\ 



2m a k [D 1 (5 ac r b + 5 ba r c ) + D 2 S cb r a } , 



(5) 



where Di and D 2 are dynamical coefficients. Their empirical values may be obtained from 
the following specific processes 



(7r p — > 7r + 7r n) — ► iT] th = 2\p2 m cr-k Di , 
(vr+p -> n + n + n) -> iT] th = V2mcr-k (D l + D 2 ) 

The pion-pole amplitude for on-shell nucleons is 



(6) 
(7) 
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cba u (p /_ p) 2_ /i 2 > 

where T^J a is the pion scattering amplitude. At tree level, it is given by 

TZa = 4 {tadtbc [{q+q'?-li 2 ]+5 bd 5 ac [{k-q'f-fi 2 ] + 6*6* [(k-q) 2 - li 2 }} 

J 7T 

and yields 

A n = _ TJWA (p'-p-k) 2 - LI 2 

f% (p'-p) 2 -Ll 2 

The contributions to A, calculated with the PV Lagrangian, are given by 



(8) 



(9) 
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the corresponding expressions for B and C are obtained by making k — > —q and A; 
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respectively, and E is 
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+ 



(12) 



(k 2 -2p'-k)(q 2 -2p-q) (q' 2 + 2p'-q')(q 2 -2p-q) (q 2 + 2p'-q)(q' 2 -2p-q') J 
For future purposes, one notes that if the PS Lagrangian were used, one would obtain 
the same structure with g A = 1 and the last line of the eq. for A would vanish. In the PS 
case, the signature of chiral symmetry are the contact interactions due to the function / in 
eq.(||), which give rise to the terms proportional to 1/m in eq. flTTD . 

In order to estimate the accuracy of the pion-production amplitude derived from eq.(0), 
we consider the contributions to the amplitudes D\ and D2, which are written in terms of 
the variables 

/i {2uj — ji) 



and 
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where uj = [/x(4m + 5/i)]/[2(m + 2/x)] and E = m + 2/i — u>. Expanding these amplitudes in 
powers of fi/m, one has 



D x = [Dl] + [D x ] 



9a 



3 a 
3 + -^ + - 

2 m 
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(17) 



D 2 = [Dl] + [D 2 ] 



^- (3 + - ±- + 

8/3 V 2 m 



(18) 



The results for the full amplitudes at threshold, namely D\ = <7a(1 + 7/x/2m)/8/ 3 and 
D 2 — — fi , yi(3 + 17/i/2m)/8/^, coincide with those obtained in the framework of chiral pertur- 
bation theory ]13[ . In order to asses the role of chiral symmetry in this problem we note that 
a Lagrangian without any symmetry, containing just a PS 7rN interaction, would give rise 
to the same pion-pole contribution and an amplitude A corresponding to just the six first 



terms in eq.(ll~D, which involve two nucleon propagators. Therefore the full elimination of 
chiral symmetry would yield D x — > g^(2)/8f^ and D 2 — > —gA(—^ + 3jt//m)/8/^, indicating 
that chiral symmetry does play a role in this problem. On the other hand, that this role 
is not as large as in the case of 7rN— > irN, where the same procedure would change one of 
the scattering lenghts by a factor of 200. Numerical results for the amplitudes are given 
in table | and show that predictions from the minimal chiral model are close to empirical 
values, although there is some room for improvement in D 2 . 



TABLE I. Subamplitudes D\ and D 2 in \x 3 units. Experimental results correspond to a best 
fit quoted in ref.[12]. 



eqs.(g§-g§) 
CHPT, eqs.(0ED 
experiment 





Dx 


Dx 


Dl 


D 2 


D 2 


1.78 


-0.91 


0.87 


-2.02 


-0.30 


-2.32 


1.81 


-0.96 


0.85 


-2.06 


-0.33 


-2.39 






0.80 






-3.20 
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In this work we are interested in the construction of the NN interaction due to the 
exchange of three pions, which is based on the amplitude T c ^ a . As indicated in fig. 3, the 
complete evaluation of this amplitude would require the calculation of a large number of di- 



agrams. However, long ago Olsson and Turner [14]] have shown that the leading contribution 



to this amplitude comes from the effective Lagrangian 

£=|4^75T^0^ 2 . (19) 



It gives rise to the following contribution to A 

A = ^ (2m- (20) 

and, as before, B and C are obtained by making k — > —q and k — > —q' respectively. This 
corresponds to the threshold amplitudes 



^-§fe(^)-f§(^)> (21 » 



- gA_ / 2m /2/A _ (2j£ 
2 8/3 V E+m Km) ~ 8fi\m> ' 



showing that the effective Lagrangian given by eq . (|19[) reproduces correctly the leading 
contribution at threshold. 



III. NUCLEON-NUCLEON INTERACTION 

The basic element in the construction of the NN potential due to the exchange of 
three uncorrelated pions (3PEP) is the corresponding Born amplitude for the process 
N(px)N(p 2 ) — ^(p'jNQog), associated with the diagrams of figure 2. Denoting this amplitude 
by F, one has 

f = u j m m aw T «° tsi ' <23) 

where the factor 1/3! is due to the symmetry of the intermediate three-pion state, A is a 
pion propagator and T^ a is the pion production amplitude for nucleon i. 
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FIG. 2. Leading contribution to the three-pion exchange potential. 
We adopt the following external kinematic variables 

W = p 1 +p 2 =p' 1 +p , 2 , (24) 

A = p[-p 1 =p 2 - p' 2 , (25) 

2=[(pi+K)-(P2+P2)]/2- (26) 
As the nucleons are assumed to be on-shell, they are constrained by 

W-z = W-A = z-A = . (27) 
For the internal variables we define 

Q = (q + q' + k)/2, (28) 

Q' = W ~ q)/2 , (29) 

so that 

k = Q - A/2 , (30) 

g = Q/2 + A/4-Q' , (31) 

q' = Q/2 + A/4 + Q' (32) 

and the condition of momentum conservation reads q + q' — k = A. 

As discussed in the previous section, the leading contribution to the amplitude T comes 
from the effective Lagrangian given by eq. (|T§D , which yields the following intermediate ef- 
fective vertex for nucleon (2) 



T!£l = i l-lp) (riyW 50 [fa® (q + q% + 8 ac r™ (q' - k), + 5 ab T^ (q - k),\ . (33) 



The corresponding expression for nucleon (1) has the same form, but is globally multiplied 
by (-1), due to the senses of flow of internal momenta. Using these results into eq.fl2"3|), we 
have 

2 



l5 uf ] V(A) 



(34) 



where 1^ is given by 



WA) = ^ / 



1 f d A Q 



(2tt)* [(Q-A/2) 2 -/i 2 ] 

d A Q' [3Q ll Q u -Q^Av/2-A fl Q„/2+27A ll A l/ /4 + AQ'^Q' V ] 
(2tt) 4 [(Q' - Q/2 - A/4) 2 - /i 2 ] [{Q'+Q/2+ A/4) 2 - ^] 

This function is evaluated in the appendix and reads 

I, v = A^A V fi 4 P(A) + 9ftt/ /i 6 I a (A) , 



(35) 



(36) 



with I p and I a given by eqs.( |A12] ) and ( |A13|) respectively. Using the nucleon equation of 
motion, one has 



F(A) 



9a 



r W. T ( 2 ) 



4m 2 /i 4 (m 75 m) (1) (m 75 m) (2) P(A) 
- (u 7 "7 6 u) (1) (u W uf ] I a (A) 



(37) 



This expression corresponds to the exchanges of pseudoscalar and axial systems. In order 
to make the strength of the interaction more transparent, we eliminate qa in favour of g, by 
means of the G-T relation, qa = gf^/m, and write 



F(A) 



/r (i). T (2)^i_j [ ( - 75u)( i) ( - 75u)(2 ) /P(A) 



/r 



Am 2 



(uYl^) {1) (u lfll5 uY z> I a (A) 



(2) 



(38) 



For future purposes, it is worth noting that the corresponding amplitude for one-pion 
exchange is 



F 7r (A) 



g 2 T W.t« 



(m7 5 m) (1) (m7 5 m) {2) 



1 



A 2 - ji 2 



(39) 
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Going to the non-relativistic limit in the center of mass frame, one has 

2 



4m 2 



/'■A<t (2) .AF(A) 



(40) 



This result allows the configuration space potential to be written as 



V 37T (x) = - — 



d 3 A 



iA-r 



4tt J (2?r) 



47T 

-/(A) 
H 



'it. tLY _/*(!). T (2) 

2m 2/2 I 4tt 



x 



(41) 



where the U(x) are integrals of Yukawa functions, written in terms of the variable x = \xr 
and given by eqs.( |A19|) and ( |A20| ). 
Using the result 

^ 2 r 7 27 60 60 



3 3 \ e 
1 + — + 



ax a 2 x 2 J x 3 



a 
~3 



W.^ 1 + JL+- + 



+ 



ax a 2 x 2 a 3 x 3 a 4 x 4 



10 45 105 105 

+ 0\2 [ 1 H 1 o — n H 5 — □ + 



ax a 2 x 2 a 3 x 3 a 4 x 4 



x" 



(42) 



with S 12 = (T^-x - o-W-crW, into egs- flAlII and (|A20| ), we obtain 



V 37T (r) 



1 ( 9H 3 V 
3 \AmPj 4tt 



L T (l). T (2) { CT (D. CT (2) + 3 (J^,)] + 0* (a;) j j (43) 



where 



l l r 1 



(4tt) 4 6 Jo Jo 



da f d~i~ia A {[3 + (l-2a) 2 ](l-2 7 ) 2 + 2[l-(l-2a) 2 ](l-2 7 ) 



7 27 60 60 \ e 
+ + 



mix) 



+ 27 - (1 2n) .} 1 • — , , , ., ; 

' \ ax a z x z a 6 x 6 arx^ 

1 1 /■! 



X 



tt.C 

3~ ' 



(44) 



(4tt) 



- [da [ d 77 a 4 f[3 + (l-2a) 2 ](l-2 7 ) 2 + 2[l-(l-2a) 2 ](l-2 7 ) 
6 jo Jo L 



. i / 10 45 105 105 \ e" a:r 
+ [27+(l-2«) 2 ] 1 + _ + _ + _ , + —-) -=-, (45) 



ax a 2 x 2 a 3 x 3 a 4 x 4 



x" 



f/ a (x) 



(4vr) 4 3 Jo Jo 

6 15 
x 1 + — + 



da / d 7 a(l — a) 7 3 (l — 7 ) a 5 



15 
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(46) 
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and 



1 



+ 



4 



V 1 -7 7[1- (1 -2a) 2 ] ' 
The structure of this result is similar to that of the OPEP, which is given by 



V*( r ) = - ( A T W. T W Ld).^) f — ] + 5 
V ; 3 V2m/ 4tt I 1 x 
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FIG. 3. Profile function for the Vq, Kf an d V a component of the three-pion exchange potential. 

The profile functions of the spin-spin and tensor components of the three-pion exchange 
potential are displayed in fig. 3, where it is possible to note that all the curves show the 
typical divergent behaviour of unregularized potentials at the origin. Therefore we assume 
that our results are realistic for internucleon distances larger than 0.7 fm, the usual bag 
radius. Inspecting the figure for the spin-spin chanel, one learns that the contribution of the 
axial component is quite small and hence the three-pion exchange potential is dominated by 
the pseudoscalar channel. In both Vss and Vr its contribution tends to add to the OPEP 
and is visible up to 1.5 fm, as shown in fig. 4. The influence of this component of the force 
over observables will be discussed elsewhere. 
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r(fm) 

FIG. 4. Profile function for the spin-spin (SS) and tensor (T) components of the (dashed 
line) and V n + V 3w (solid line). 



In order to produce a feeling for the structure of the functions U(x), in the appendix 
we have evaluated approximately the integrals in eqs.(^4|446]) and obtained the following 
asymptotic results (x oo) 



U$(x) 
Ul{x) 
U a {x) 



7T 80 A 3 13 10 10 \ e 



(4tt) 4 y/3 V x 3x 2 3x 3 9x 4 / x 4 



7T 80 



80 / 4 20 16 16 \ e 



(4tt) 4 ^3 V x 3x 2 3x 3 9x 4 / x 4 



-3.c 



7T 16 / 2 5 5 \ e 

1 + - + 7T^ + ^ 



(4tt) 4 3^ V x 3x 2 9x 3 J x 5 



-3x 



(49) 
(50) 
(51) 



which are compared with the exact ones in fig. 5. 
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r(fm) 

FIG. 5. Ratios between the approximate expressions (49-51) by the corresponding functions 
(44-4^) - all ratios are identical. 

A last point we would like to address here concerns the nature of the force in the chiral 
limit. The potential given in eqs.(ff3|-|4"TD incorporates two kinds of approximations. The first 
of them is associated with the assumption that C, eq.(|H?D, represents the leading contribution 
to the NN7T7T7r vertex. The other one is related to the non-relativistic limit taken in eq. ([H]) . 
On the other hand, no approximations besides the neglect of contact interactions were 
performed in the calculation of the three-pion propagator represented by the functions /(A). 
Therefore the corresponding configuration space expressions, given by eqs.([H]-46) also do 
not contain approximations and can be used to evaluate the form of the interaction in the 
chiral limit. The strength of V 3lT (r), as given by eq . (|43|) , is proportional to /i 7 . Recalling 
that x = fir, we obtain the following results when \i — > 0: h 7 Uq(x) — > 140/[(47r) 4 r 7 ], 
fj, 7 U2{x) — > 245/[(47r) 4 r 7 ] and fi 7 U a (x) — > — 5/[(47r) 4 r 7 ]. Thus, the three-pion exchange 
NN potential survives in the chiral limit, where it has the form 

N " i ,1 

- . (52) 



V^(r) 



.(i). -(2) 



3 \4:mfl J (4tt) e 



125 <t (1) -<t (2) + 245 S 
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APPENDIX A: INTEGRALS 

In this appendix we evaluate the integral given by eq. (|35|) , using the following results 

d A Q 1 



da 



(2tt) 4 [(Q-K/2f-^) [{Q+K/2)*-tf?\ 
Po-lnf l--^2 



p 2 £ 2 



(4tt) 2 Jo 

A " SA:/ ''° I (2 [(()+ K/2) 2 — fi 2 t, 2] 



(Al) 



d*Q 



Ql 



(4tt 



,2 A* 



A„ / da 



1 - 2a 



Po-ln 1- 



A 2 



(A2) 



Q iiQv 



d 4 Q 

(2tt) 4 [(Q-A/2) 2 -p 2 ] [(Q + K/2) 2 -^e 



(4tt) 



K^K V da 



K /l-2a 



P °- ln ( 1 -^ 



p 2 ri S 2 [l-(l-2a) 2 ] / A 



4 



p 2 £ 2 



pi-ln 1- 



A 2 ' 



(A3) 



where 



l2 4[a + (1 - a)i 2 \ 



(A4) 



[1 - (1 - 2a) 2 ] 

Po and pi are constants associated with the dimensional regularization procedure. In order 
to perform the integrations, it is convenient to use the following representation for the 
logarithm 

A 2 N 



In 1 



p 2 £ 2 



L d/3 G + K*-J>T,yp) ■ 



(A5) 



Quite generally, constants appearing in these results correspond to contact interactions, 
since they do not depend on A. As we are interested in the long range part of the potential, 
these constants will be neglected in the sequence and we write 
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(A6) 
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The integral in Q' can be performed using these results and one has 

d 4 Q' [3Q li Q I/ -Q ll A v /2-A IJi Q v /2+27A li A v /A + AQ'^Q' U ] 



(2tt)4 [(g'_g/ 2 -A/4) 2 - /U 2 ][(g' + g/2 + A/4) 2 - /U 2 ] 
(3g M g,-g^A,/2- A^/2+27A At A,/4)/ / u 2 X(g+A/2; //, 1) 



+ 4I F (Q+A/2;/i,l) 

= - (47F r^r^ ^ [(g+A/^-^A 2 ] {( 3+ ( i - 2a ) 2 ] 

- [l-(l-2a) 2 ](g M A ! , + A M g,)/2+ [27+(l-2a) 2 ] A M A,/4 

- A 2 [l-(l-2«) 2 ](l-/3)/i 2 ^/2} , 



(A9) 



where 



A 2 = 



/3[l-(l-2a) 2 ] • 



(A10) 



The function is then given by 

/ (A ) = i f*9- I I' 

^ y ' 3! J (2tt) 4 \{Q- A/2) 2 -/z 2 ] m " 

i l , r 1 . r l A 2 



-/x 2 | o rfa | Q rf/3-{[3+(l-2«) 2 ]X^(A,/i,A) 



(4tt) 2 3! ^ Jo Jo '~ (3 

- [1 - (1 -2a) 2 ] [A„X„(A, ii, A) + A„X„(A, ^ A)]/2 + [27+ (1 -2a) 2 ]X(A, ^, A) A M A,/4 

- A 2 [l-(l-2a) 2 ](l-/3) X(A >A1 , A) /i 2 ^/2} 

= A M A, /i 4 F(A) + ^ /i 6 /"(A) , (All) 
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where 



+ 2[l-(l-2 a ) 2 ](l-2 7 ) + |27+(l-2 a ) 2 ]} [A2 ^ 2(p] , (A12) 
+ f> 2 [3+(l-2 a ) 2 ][l-(l-2i) 2 ](l- e )} |A2 _^ 2g2] , (A13) 



with 



4[ 7 + (l- 7 )A 2 ] 



^ = 4^ 44 . (A14) 

c[l- (1-27)*] 1 J 

Results presented in this appendix are covariant. On the other hand, the calculation of 
the potential is performed in the centre of mass of the NN system and we use A = (0; A) — > 
A 2 = -A 2 . 

The potential in configuration space is determined by the functions U(x), given by 

^7/S e - A - r/(A) - (A15) 

where x = fir. Using the result 

J (2tt) 3 A 2 + /i 2 # 2 4tt x ' { ' 

we have 

-Ox 

+ 2[l-(l-2a) 2 ](l-2 7 ) + [27+(l-2a) 2 ]} , (A17) 

+ # 2 [3 + (l-2«) 2 ][l-(l-2 7 ) 2 ](l- e )| — . (A18) 



X 

The integrations in e and j3 can be performed analytically and we have 
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U p (x) 



l l r 1 



da f d-f-fa 2 {[3 + (l-2a) 2 ](l-2 7 ) 2 + 2[l-(l-2a) 2 ](l-2 7 ) 



(4tt) 4 6 Jo Jo 
[27+(l-2«) 2 ]}fl + — + 3 

> \ n T. 



ax a?x 2 J x 3 



(A19) 



U a (x) = - - J- | ['da fd 1 a(l - a) 7 3 (l - 7) a 5 

(47T) 4 6 JO JO 



x 1 



6 15 
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ax a?x 2 a 3 x 3 



x' 



(A20) 



where 



+ 



(A21) 



V 1 -7 7[1- (1 -2«) 2 ] ' 
The integrals over a and 7 can be evaluated approximately for large values of x. In this 
case, the exponential has a sharp minimum for a = 1/2,7 = 2/3 and varies very rapidly 
around it. Thus all the elements in the integrand but the exponential may be taken as 
constants and we have 



(4tt) 

U a (x) = - t— - 16 ( 1 ~ 



(4tt) / 



x 



5 5 
3x 2 9x 3 



1 r l r 1 

da d^f e 



x^ Jo 



(A22) 
(A23) 



In order to perform the last integral we first use a new variable u, related to a by 

1/2 



1 

a — - 
2 



1 



1 



_4 4 + 2 7M 2 v /^ + 7 ^ 
and then another variable v, related to 7 by 



(A24) 



7 = 



3 + (3 + v 2 f =f V48t; 2 + 44w 2 + 12v 6 + 



(A25) 



2(3 + t; 2 ) z 

where the (— ) and (+) signs refer to the intervals < 7 < 2/3 and 2/3 < 7 < 1 respectively. 
We then obtain 

-3x 



I da d*y 
Jo 



7r e 
3^ x 



(A26) 
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